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Abstract. We develop a general method to bound the spreading of an entire 
wavepacket under Schrodinger dynamics from above. This method derives 
upper bounds on time-averaged moments of the position operator from lower 
bounds on norms of transfer matrices at complex energies. 

This general result is applied to the Fibonacci operator. We find that at 
sufficiently large coupling, all transport exponents take values strictly between 
zero and one. This is the first rigorous result on anomalous transport. 

For quasi-periodic potentials associated with trigonometric polynomials, we 
prove that all lower transport exponents and, under a weak assumption on the 
frequency, all upper transport exponents vanish for all phases if the Lyapunov 
exponent is uniformly bounded away from zero. By a well-known result of 
Herman, this assumption always holds at sufficiently large coupling. For the 
particular case of the almost Mathieu operator, our result applies for coupling 
greater than two. 



1. Introduction 

If H is a self-adjoint operator on a separable Hilbert space Tt , the time-dependent 
Schrodinger equation of quantum mechanics, idtip = Hip, leads to a unitary dy- 
namical evolution in Ti., 

(i) m = e - itH m- 

Of course, the case of main interest is where Ti. is given by L 2 (K d ) or £ 2 (Z d ), H is 
a Schrodinger operator of the form —A + V, and ip(0) is a localized wavepacket. 

Under the time evolution Q), the wavepacket ip(t) will in general spread out 
with time. Finding quantitative bounds for this spreading in concrete cases and 
developing methods for proving such bounds has been the objective of a great 
number of papers in recent years. 

These questions are particularly complicated and interesting in situations where 
singular continuous spectra occur. The discovery of the occurrence of such spectra 
for many models of physical relevance has triggered some activity on the quantum 
dynamical side. For models with absolutely continuous spectra, on the other hand, 
quantum dynamical questions have been investigated much earlier in the context of 
scattering theory and have been answered to a great degree of satisfaction. Models 
with point spectrum or, more restrictively, models exhibiting Anderson localiza- 
tion, have also been investigated from a quantum dynamical point of view. While 
it is not true in general that spectral localization implies dynamical localization 
|2fil IT3| , it has been shown that the three main approaches to spectral localiza- 
tion, the fractional moment method 0EJO3, multi-scale analysis |^H1 H3 UM > an d 
the Bourgain-Goldstein-Schlag method |121 137| . yield dynamical localization as a 

byproduct |2 E2 E2 EU EU . 
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Quantum dynamics for Schrodinger operators with singular continuous spectra 
is much less understood. The two primary examples motivated by physics are 
the Harper model, also known as the almost Mathieu operator at critical cou- 
pling, which describes a Bloch electron in a constant magnetic field with irrational 
flux through a unit cell, and the Fibonacci model, which describes a standard 
one-dimensional quasicrystal. Both models have been studied heavily in both the 
physics and mathematics communities. 

The standard approach to dynamical questions is via the spectral theorem. Re- 
call that each initial vector ip(0) = tp has a spectral measure, defined as the unique 
Borel measure obeying 



for every measurable function /. Here, (•, •) denotes the scalar product in Ti. A very 
important discovery of Guarneri |89II4()| . which has been extended by other authors 
[T4l rm H2l EH) 7 was that suitable continuity properties of the spectral measure 
dfi^ imply lower bounds on the spreading of the wavepacket. Such continuity 
properties follow from upper bounds on the measure of intervals, n^([E — e, E + 
e]), E € c(-ff), £ — * 0. Later on, extensions and refinements of Guarneri's method 
were developed |El E2l EH- which allowed the authors of those papers to obtain 
more general dynamical lower bounds that take into account the whole statistics 
of H^([E — e, E + e]), E € K. In particular, lower bounds on ^,([E — e, E + e]) 
for E from some set of positive Lebesgue measure imply dynamical lower bounds 
in the case of both singular continuous and pure point spectrum . Better lower 
bounds can be obtained if one has information about both the measure of intervals 
and the growth of the generalized eigenfunctions u^(n, E) |58l I71| . 

In the case of Schrodinger operators in one space dimension, the required infor- 
mation on the spectral measure (and on generalized eigenfunctions) is intimately 
related to properties of solutions to the differential/difference equation Hu = Eu 
PUl l2"51 0^ E2 E|. An analysis of these solution has been carried out in a 
number of concrete cases ^Hl EEH 023 EU E2l EU , giving rise to explicit lower bounds 
on spreading rates for (generalizations of) the Fibonacci model, sparse potentials, 
and (random) decaying potentials. Thus, in one dimension, there is a three-step 
procedure for obtaining lower bounds on spreading rates for a given model: prop- 
erties of the model ~> solution estimates ~> bounds on the measure of intervals 
H^([E — e, E + £•]) ~~> dynamical lower bounds. 

There is a second approach to dynamical lower bounds in one dimension, which 
is based on the Parseval formula, 1 



(2) 2tt / e- 2t / T \(e- ltH 8 1 ,8 n )\ 2 dt= \((H - E - i)" 1 ^, S n )\ 2 dE. 



This method was developed in [231 1241 . It has the advantage that it gives dy- 
namical bounds rather directly without any knowledge of the properties of spectral 
measure (although additional information on spectral measure allows one to im- 
prove the results [211 1121 )• The input required for this method is based on upper 
bounds for solutions corresponding to some set of energies, which can be very small 
(non-empty is sufficient 123 021 )• This method is the basis for the results in |21H53| . 

lr The formula was used earlier by Kiselev, Killip, Last in their study of dynamical upper 
bounds on the slow part of the wavepacket 1571 . which will be discussed momentarily. 
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Moreover, a combination of this approach with the traditional one (based on the 
spectral theorem) leads to optimal dynamical bounds for growing sparse potentials 
|72| and the best dynamical lower bounds for the Fibonacci operator known to this 
date [25] . 

To a certain extent, there is a fairly good understanding of how to prove dy- 
namical lower bounds, particularly in one space dimension, as a consequence of the 
results described above. The problem of proving upper bounds, on the other hand, 
is much harder and has essentially eluded the researchers in this field up to this 
point. The known lower bounds for models with singular continuous spectra do not 
exclude ballistic motion and hence are not able to distinguish them from models 
with absolutely continuous spectra. However, ballistic motion is not expected to 
occur for the two main models of interest, the Harper model and the Fibonacci 
model. While we explain the terminology below, let us mention at this point that 
the motion is expected to be almost-diffusive in the Harper model, while it is ex- 
pected to be anomalous (i.e., neither localized, nor diffusive, nor ballistic) in the 
Fibonacci model. There is a lot of numerical and heuristic evidence for these claims 
in the physics literature. To prove such a statement, one needs two-sided estimates 
for the dynamical quantities associated with these models. 

Proving upper bounds is so hard because one needs to control the entire 
wavepacket. In fact, the dynamical lower bounds that are typically established 
only bound some (fast) part of the wavepacket from below and this is sufficient for 
the desired growth of the standard dynamical quantities. It is of course much easier 
to prove upper bounds only for a portion of the wavepacket and Killip, Kiselev and 
Last succeeded in doing so jHZ]- Namely, their work provides explicit criteria for 
upper bounds on the slow part of the wavepacket in terms of lower bounds on so- 
lutions. They applied their general method to the Fibonacci operator. Their result 
further supports the conjecture that this model exhibits anomalous transport. 

It is the purpose of this paper to develop a general method connecting solution 
properties to dynamical upper bounds and to use this method to prove dynamical 
upper bounds for the Fibonacci operator and the almost Mathieu operator. 

The general method is based on the Parseval formula (J5J, which has already 
proved its utility when studying dynamics. We are able to bound the entire 
wavepacket from above, provided that suitable lower bounds for solution (or rather 
transfer matrix) growth at complex energies are available. 

We apply this method to the Fibonacci operator and establish anomalous trans- 
port at sufficiently large coupling. In combination with the known lower bounds, 
our result shows that the large coupling asymptotics of the dynamics follows a law 
predicted by Abe and Hiramoto in 1987. 

Finally, we prove upper bounds for the almost Mathieu operator and extensions 
thereof at sufficiently large coupling. It is shown that positive Lyapunov exponents 
imply a weak form of dynamical localization in a sense that will be made precise 
below. This result complements a stunning observation by Last who proved al- 
most ballistic transport on a sequence of time scales for the same model in the 
case of some Liouville frequencies. Our result shows that on other sequences of 
time scales, one has essentially no transport. Thus, the almost Mathieu operator 
with a Liouville frequency exhibits quite interesting dynamics. We also show that 
the phenomenon discovered by Last is impossible if the frequency satisfies a weak 
Brjuno-type condition. 



4 



DAVID DAMANIK AND SERGUEI TCHEREMCHANTSEV 



Acknowledgments. It is a pleasure to thank Nikolai Makarov and Mihai Stoiciu for 
useful conversations. 

2. Description of the Results and Overview 
Consider a discrete one-dimensional Schrodinger operator, 

(3) [Hip](n) = if>(n + 1) + ^{n - 1) + V(n)^(n), 

on £ 2 (Z) or ^ 2 (Z+), where Z+ = {1,2, . . .}. In the case of l 2 {1+), we will work 
with a Dirichlet boundary condition, ip(0) = 0, but our results easily extend to all 
other self-adjoint boundary conditions. 

A number of recent papers (e.g., gSl EH HS1 EH E2 E31 EZHZH E3 ) were devoted 
to proving lower bounds on the spreading of an initially localized wavepacket, say 
ip = 5\, under the dynamics governed by H, typically in situations where the 
spectral measure of <5i with respect to H is purely singular and sometimes even 
pure point. 

A standard quantity that is considered to measure the spreading of the 
wavepacket is the following: For p > 0, define 

(4) (|X|? i )(T)=^|n|Mn,T), 

n 

where 

9 r°° 

(5) a(n,T) = ~ J e- 2t / T \(e~ itH S 1 ,S n )\ 2 dt. 

Clearly, the faster (\X\g )(T) grows, the faster e~ ltH 8\ spreads out, at least averaged 
in time. One typically wants to prove power-law bounds on )(T) and hence 

it is natural to define the following quantity: For p > 0, define the lower transport 
exponent /3J" (p) by 

and the upper transport exponent fjg (p) by 

log(\Xf Si )(T) 
^»=hmsup plQgT . 

Both functions /3 S (p) are nondecreasing. 

Another way to describe the spreading of the wave-function is in terms of prob- 
abilities. We define time-averaged outside probabilities by 

P(N,T)= Yl a ^ T )- 

\n\>N 

Following [5J|, for any a S [0, +oo] define 2 
6 S (a) = -lrminf — 

T^oo log 1 

and 

(7) sn^-^up ^f'" 2 ^ 

T^oo log T 



2 We take T a - 2 so that P(T° - 2, T) = 1 for all T. 
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For every a, < S + (a) < S~(a) < oo. 

These numbers control the power decaying tails of the wavepacket. In particular, 
the following critical exponents are of interest: 

(8) af = sup{a > : S ± (a) = 0}, 

(9) af =sup{a > : S ± (a) < oo}. 

We have that < aj < a~ < 1, < af < af < 1, and also that aj < af, 
a~ < af. One can interpret af as the (lower and upper) rates of propagation 
of the essential part of the wavepacket, and af as the rates of propagation of the 
fastest (polynomially small) part of the wavepacket; compare [33] . In particular, 
if a > af, then P(T a ,T) goes to faster than any inverse power of T. Since a 
ballistic upper bound holds in our case (for any potential V), Theorem 4.1 in j.'i.'ij 
yields 

lim 8f i (p)=af 

and 

lim Pf 1 (P) = <- 
In particular, since (p) are nondecreasing, we have that 

(10) 13% {jp) < af for every p > 0. 

When one wants to bound all these dynamical quantities for specific models, it is 
useful to connect them to the qualitative behavior of the solutions of the difference 
equation 

(11) u(n + 1) + u(n - 1) + V(n)u(n) = zu(n) 

since there are effective methods for studying these solutions. Presently, the 
known general results are limited to one-sided estimates of the transport expo- 
nents. Namely, as already alluded to in the introduction, a number of approaches 
to lower bounds on /3 S (p) have been found in recent years. 

It should be stressed that there were no general methods known to bound af , a f , 
or (p) non-trivially from above. In the present paper we propose a first general 
approach to proving upper bounds on af (which in turn bound af and ^{p) for 
all p > from above, as well). This approach relates the dynamical quantities 
introduced above to the behavior of the solutions to the difference equation l)ll|l 
for complex energies z. To state this result, let us recall the reformulation of 
(|llfl in terms of transfer matrices. These matrices are uniquely determined by the 
requirement that 

U{n+1) )-<P(n z)( 
u{n) ) ~ q?[n ' Z > { u(0) 

for every solution u of (|11|) . Consequently, 

'T(n,z)---T(l,z) n>l, 

(12) <f>(n,z) = {Id n = 0, 

M)]- 1 ."^,*)]- 1 n<-l, 

z — V(m) —1 
1 




where 

T(m, z) 
We have the following result: 
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Theorem 1. Suppose H is given by lj2J), where V is a bounded real-valued function, 
and K > 4 is such that a (H) C [— K+ 1, K — 1] . Suppose that, for some C € (0, oo) 
and a € (0, 1), we have 

(13) |^ (i<^gr«. Il $ ( n ' B + t) f) = °( T_m ) 
and 

(14) ^ma^J^^E+^ll 2 ) dE = 0(T- m ) 

for every m > 1. TTien 

(15) a+ < a. 
In particular, 

(16) P + (p) < /or every p > 0. 

Remarks, (a) If the conditions of the theorem are fulfilled for some sequence of 
times Tfe — > oo, we get an upper bound for a~. 

(b) The statement of the theorem follows from upper bounds for outside probabil- 
ities described in Theorem below. 

The next issue we want to address is the stability of the transport exponents 
with respect to suitable perturbations of the potential. For example, it is reasonable 
to expect that they are invariant with respect to finitely supported perturbations. 
However, such a property has not been established yet. Moreover, the approaches 
to lower bounds for transport exponents that are based on dimensional properties 
of spectral measures are not suited to prove this invariance. It is well known that 
these dimensions are rather sensitive with respect to finite rank perturbations. As 
pointed out in [23 E], the approach developed in those papers is stable in the sense 
that if it can be applied to a certain model, it can also be applied (and yields the 
same bounds) to finitely supported (or even suitable power-decaying) perturbations. 

Of course, it is even more desirable to prove that the transport exponents (as 
opposed to one-sided bounds) are invariants. The following result establishes this 
fact: 

Theorem 2. Let Hi, H2 be two operators of the form (PJ with bounded potentials 
and let K > 4 be such that a (Hi. 2) C [-K + 1,K — 1]. Denote the corresponding 
transfer matrices by $1, $2 and the corresponding transport exponents by ctf u , ctj u . 
Assume there exists A > such that for all E £ [—K,K], < e < 1, |n| < 1/e, 

(17) e A \\$i(n,E + ie)\\ < \\® 2 (n, E + ie)\\ < e" A ||$i(n, E + ie)\\. 
Then, af >u = af u . 

The bounds (|17J) clearly hold if the potentials differ only on a finite set. Thus, 
Theorem |21 yields the invariance mentioned above. For example, assume that dy- 
namical localization holds for the operator Hi. That means, in particular, that 
@t 81 (p) = f° r ever y P > an d thus a~l_ u = 0. Theorem yields a^ u = 0, so that 
P2 Si (p) = — some weak form of dynamical localization for H 2 . This result should 
be compared with that of |26j about semi-stability of dynamical localization. 
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In most cases where nontrivial dynamical lower bounds have been proven, the 
results obtained imply af > and a f = 1. This is the case for sparse potentials 
[331 I72| . random decaying potentials the Thue-Morse Hamiltonian and 
random polymer models |53| . Although it has not yet been proven, it is expected 
that af < 1 for some of these models (for sparse potentials, however, it was shown 
in [E] that af = 1). 

On the other hand, there are models for which it is expected that a~ may take 
values other than zero (this is the case when one has dynamical localization) or one 
(ballistic transport at least for p> 1). Such a situation is called anomalous trans- 
port by some authors. More restrictively, a f = 1/2 is sometimes called diffusive 
transport and anomalous transport refers to a situation where one does not have 
ballistic transport, diffusive transport, or dynamical localization. 

Two very prominent examples are given by the Fibonacci operator and the almost 
Mathieu operator. We will apply our general result, Theorem^ to both of them. 

The Fibonacci operator has been studied for decades by mathematicians and 
physicists. The potential is given by 

(18) V{n) = X X [i-e,i){n6 modi), 8 = ^—^-. 

This potential belongs to the more general class of Sturmian potentials, given by 

(19) V{n) = \x[i-8A){nO + w modi) 

with general irrational 9 £ (0, 1) and arbitrary uj E [0, 1). These sequences provide 
standard models for one-dimensional quasicrystals. 3 Early numerical and heuris- 
tic studies of the spectral properties of the Fibonacci model were performed by 
Kohmoto, Kadanoff, Tang [^j and Ostlund, Pandit, Rand, Schellnhuber, Siggia 
|K4j . It was suggested that the spectrum is always purely singular continuous. This 
was rigorously established by Siito for the Fibonacci case [SO [70] and by Bellissard, 
Iochum, Scoppola, Testard JT] and Damanik, Killip, Lenz in the general Stur- 
mian case. Abe and Hiramoto studied the transport exponents for the Fibonacci 
model numerically [Tl 144) . They found that they are decreasing in A and behave 
like 

(20) af, af - ' 



log A 

as A — > oo. (Here and in the following, we write / x g if C~ x g < f < Cg for 
some C > 1.) See also [3PJ|SS] for other numerical studies of Fibonacci quantum 
dynamics. 

The general approaches to lower bounds for the transport exponent described 
above have all been applied to the Fibonacci Hamiltonian (and some Sturmian 
models). The best lower bound for af was obtained by Killip, Kiselev, and Last 
in [571 ■ It reads 4 

2k 



> 



C(A) +K + 1/2' 



^See Ifi (il for the discovery of quasicrystals and |H| for surveys on the mathematical theory of 
quasicrystals in general and the role of the Fibonacci operator in particular. 

4 The expression for ((\) is given in terms of its large A asymptotics which is of main interest 
here. See [24 57 for its values at small A. 
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where 

0.0126 



and C(A), chosen so that one can prove a result like 

J2\Mn,E)\\ 2 <CL 2 ^ +1 

n=l 

for energies in the spectrum of H (our definition differs from that of [FT]), obeys 

31ogvo , 
C(A) = ^=-y (log A + O(l)) . 

log (^tl) 

This shows in particular that af admits a lower bound of the type (|20J) . 
The best lower bound for a~ was found in |24| . where it was shown that 

1 



C(A) + 1 

In terms of the exponents (3J (p), the best known lower bounds are (see |25j) 

B-(n) > J (p+i)(C P (A)+ k +i/2) P< 2 C(A) + 1, 
^"h*! P>2C(A) + 1. 

We also want to mention work on upper bounds for the slow part of the 
wavepacket by Killip, Kiselev, and Last, |S7]. More precisely, they showed that 
there exists a S £ (0, 1) such that for A large enough (so that p(X) defined in 121|) 
below is less than one), 

P(C 2 T pW ,T) <l-6. 

Here, 

6 log ^+1 

(21) p(A) = ° g ... 2 . 

logC(A) 

and 



(22) w- '-^t^ 

See j^| Theorem 1.6. (i)] . However, this result does not say anything about the fast 
part of the wavepacket, and in particular, no statement for any of the transport 
exponents can be derived. 

With the help of Theorem^we can prove upper bounds for a+ for the Fibonacci 
model at sufficiently large coupling. These upper bounds show that l|2rj|) is indeed 
true. 

The precise result is as follows: 

Theorem 3. Consider the Fibonacci Hamiltonian, that is, the operator J3J| with 
potential (|18fl . Assume that A > 8 and let 
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with £(A) as in J22J. Then, 

(23) a+ < o(A), 
and hence 

(24) /3 + (p) < a(A) /or euery p > 0. 

One can observe that a(X) < p(X), where p(X) is the number given in (j21(l — the 
power for which [HZj proved upper bounds for the slow part. Note that £(A) = 
A + 0(1) as A — > co. Moreover, 

a(8) = — , 2 w 0.876 
log 3 

and a(A) is a decreasing function of A for 5 A > 8. Thus, we establish anomalous 
transport for the Fibonacci Hamiltonian with coupling A > 8 and confirm the 
asymptotic dependence of the transport exponents on the coupling constant A 
that was predicted by Abe and Hiramoto. We emphasize again that this is the first 
model for which anomalous transport, in the sense that < aj < a+ < 1, can be 
shown rigorously. 

The almost Mathieu operator is given by J3J with potential 

(25) V(n) = A cos(2n{n9 + w)) 

with 9 irrational and u> E [0, 1). The almost Mathieu operator in the critical case 
A = 2 is also called the Harper model, 

(26) V{n) = 2 cos(27r(n6» + lo)) 

Operators with potential i|25|) have been studied extensively by mathematicians 
and physicists for decades. The spectral theory is almost completely understood. 
For almost every frequency 9 and phase u>, one has purely absolutely continuous 
spectrum for < A < 2, purely singular continuous spectrum for A = 2, and pure 
point spectrum with exponentially decaying eigenfunctions (Anderson localization) 
for A > 2. 

On the other hand, for 9's that are very well approximated by rational numbers, 
one can use Gordon's Lemma j2U to prove absence of eigenvalues so that localization 
fails for these 9's and all w's when A > 2. But even for a Diophantine 9, absence 
of eigenvalues holds generically, so that for a dense Gg set of w's, one has again no 
eigenvalues |54|. To a certain extent, the situation was mollified by Jitomirskaya 
and Last |221, who proved that for A > 2, all irrational frequencies and all phases, 
the spectral measures are zero-dimensional. 6 

On the dynamical side, the result of Last that for < A < 2, every irrational 9, 
and every 7 uj, there exists some absolutely continuous spectrum |61) . implies ballistic 
motion in this coupling regime. In the regime A > 2, it is of course desirable to 
establish dynamical localization whenever spectral localization holds. For results in 
this direction, see |30l !32llB17| . As a result, dynamical localization holds for A > 2, 
Diophantine 9, and almost every to. 



5 The minimum value, Ao, can be chosen slightly smaller than 8. By monotonicity, all we need 
to require from Ao is that £(Ao) > 1 and o(Ao) < 1. 

6 This was proved earlier in the special case of Liouville frequencies by Last 1621 . 

7 Last proves this for almost every ui; it holds for every w by a result of Last and Simon |63l . 
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On the other hand, it was not clear what should be expected for the dynamical 
quantities (p) = /3^~ (p, A, 9, u) in cases of exceptional frequencies and/or phases 
at coupling A > 2. Since the Guarneri-Combes-Last bound is one-sided, zero- 
dimensionality of spectral measures has no dynamical implications. In fact, Last 
proved the following result in /3^(2, A, 6, iS) = 1 (and hence a+(A, 6, (S) = 1) 
for every A, uj, and a suitable (A, w)-dependent Liouville frequency 9. Del Rio, 
Jitomirskaya, Last, and Simon constructed 9 such that (p, 3, 9, lo) = 1 for every 
p > 0,lo G [0, 1) [2U (their proof given for p = 2 can be easily generalized to any 
p > 0). Thus, for such a value of 9, a/~(3, 9, lo) = 1. Later on, it was shown in (231 
that given A > 2, one has a/" (A, 9, lo) = 1 for all 9's from some dense Gg set. 

Thus, in these exceptional cases, one can only hope to bound a~ (A, 9, lo) from 
above non-trivially. 

Recall the continued fraction expansion of 9, 



a 



ai 



1 

a 2 



a 3 H 

with uniquely determined integers ao £ Z and a n > 1, n > 1; see |5(j) for background 
information. Truncation of this expansion gives rise to best rational approximants 
Pk/qk of 6, 

Pk 1 
— = ao 



qk 1 
a\ H 



0,2 



1 

(>k 



The following condition on 9 will be useful below: 

(27) lim l0ggfc+1 = 0. 

fe^oo qk 

This condition is satisfied for Lebesgue almost every 9; see, for example, [SBJ Theo- 
rem 31]. It is related to, but much weaker than, the Brjuno condition, which states 
that (logqk+i)/qk € t • The Brjuno condition is crucial, for example, in Yoccoz' 
work on analytic linearization of circle diffeomorphisms [ZZj- Every Diophantine 
number satisfies the Brjuno condition and the inclusion is strict. 

We have the following dynamical results for the almost Mathieu operator at 
super-critical coupling: 

Theorem 4. For the almost Mathieu operator with coupling A > 2, any irrational 
9, and any uj e [0,1), we have a~(\,9,u>) — (and hence /3^~ (p, A, 9, lo) = for 
every p > 0). If in addition 9 obeys l|27|l . then (A, 9,u>) — for every lo S [0, 1) 
(and hence (p, A, 6, u>) = for every p > 0) . 

Remarks, (i) In particular, for any given coupling A > 2, any phase lo, and an 
associated Liouville frequency 9, there is a sequence of time scales, for which one 
has almost ballistic transport, but there is essentially no transport along certain 
other sequences of time scales. Thus, the Liouville case may exhibit interesting 
dynamical phenomena. 

(ii) Under the rather weak condition (|27|l we get a certain form of dynamical local- 
ization. While this result is slightly weaker than the one established by Germinct 
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and Jitomirskaya under a Diophantine assumption on 9 [22], it has several advan- 
tages: 

• The assumption (|27[) is much weaker than the one from |32]; compare | 56l 
Theorems 9 and 13]. As mentioned above, any standard Diophantine con- 
dition implies the Brjuno condition, which in turn implies <|27|) . and both 
inclusions are strict. 

• Our result holds for all phases uj, while Germinet and Jitomirskaya consider 
a phase average and, in particular, lose control over zero measure sets. As 
mentioned above, even under a Diophantine condition, there are exceptional 
values of uj for which one has no eigenvalues. Germinet and Jitomirskaya 
cannot make any statement for these ui's, while we can handle every uj. 
In a way, our result bears some similarity with a dynamical upper bound 
established by del Rio, Jitomirskaya, Last, Simon in |261 Theorem 8.1] 
who considered exceptional couplings in a rank-one perturbation situation: 
While there are exceptional phases, for which spectral localization fails, the 
failure of (the strong form of) dynamical localization cannot be too severe. 

• Our proof is much shorter. It only relies on positive Lyapunov exponents 
and does not need to establish spectral localization first. Note that pT2] 
relies on |121 1371 148| . while our proof is self-contained up to a proof of 
positive Lyapunov exponents, which can be proved in a few lines 02]; see 
also [E3 pp. 199-200]. 

We will prove Theorem 0] in a more general form. Let us consider an operator 
H = Hg tU , of the form © whose potential is given by 



where / is a trigonometric polynomial, 9 is irrational, and u> G [0, 1). 

The Lyapunov exponent, whose existence follows from Kingman's subadditive 
ergodic theorem, is given by 



(28) je{z)= Mm - / log\\$ n {z,6,w)\\du> = inf - / ]og\\$ n (z,8,u)\\du. 



One gets the same quantity if one considers n < — 1. 

Theorem 5. Let Hg jLU be as just described. Assume that there is a number L > 
such that je(z) > T for every z£C. Then a~{9,uj) = for every uj £ [0, 1). If in 
addition 9 obeys 127|) . then a^{9,uj) — for every uj G [0, 1). 

A classical result of Herman 02] states that if f(uj) — Acos(27rw), then 



Thus, Theorem implies Theorem 0] 

For a general non-constant trigonometric polynomial /, Herman's argument 
shows that the assumption of Theorem [5] holds for sufficiently large leading co- 
efficient. Thus, we may state the following 

Theorem 6. Suppose g is a non-constant trigonometric polynomial and f = Xg. 
Then, there exists \q such that for every X > Xq, every irrational 9, and every 
uj € [0,1), we have a~{X,9,u>) = 0. Moreover, if 9 obeys Q27I1 . then o£(\,6,uj) = 
for every X > Xq and uj 6 [0, 1). 



V(n) = Vg, u {n) = f(n9 + uj), 





le{z) > log — for every 9 e K \ Q, z G C. 
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Note that Theorems 0] and El combined with the Guarneri-Combes-Last bound, 
provide a new proof of the zero-dimensionality of spectral measures since we do not 
use the Jitomirskaya-Last inequality from |51| on which the proof in |52| was based. 

We believe that one can extend Theorem and hence Theorem HO to the case 
of / = Xg with g analytic and A sufficiently large. Positivity of the Lyapunov 
exponent was shown in this case by Sorets and Spencer [HEJ; see also ^21 EH]- Zero- 
dimensionality of spectral measures for analytic potentials and large coupling is a 
result of Jitomirskaya and Landrigan |49| . 

Given Germinet- Jitomirskaya [32], Last [S2], and Theorem 0] above, we have a 
good understanding of the dynamical picture in the super-critical case, A > 2. Due 
to the existence of absolutely continuous spectrum in the sub-critical case 61 , 
< A < 2, this leaves the intriguing problem of studying the critical case, A = 2, 
that is, when the potential is given by l|26|) . There are no rigorous results to this 
date. Based on numerics |221 ESI EH1 EH| j h is expected that one has almost diffusive 
transport in the sense that the transport exponents are close, but not equal, to 1/2. 

The organization is as follows. In Section 3 we prove upper bounds for outside 
probabilities in terms of the norms of transfer matrices at complex energies; see 
Theorem These bounds immediately imply Theorem ^ a bove. We also express 
the numbers in terms of transfer matrix norms (see Theorem |SJ and prove 
Theorem [21 In Section 4 we study the Fibonacci Hamiltonian. Suitable bounds 
on transfer matrix norms are obtained by studying the traces of these matrices. 
Since they obey recursive relations, we have to investigate a dynamical system, the 
so-called trace map. As opposed to earlier studies of the trace map, we need to 
regard this as a complex dynamical system and use ideas from conformal mapping 
theory, the Koebe Distortion Theorem in particular, to prove the required estimates. 
Once these estimates have been established, Theorem [JJ follows from Theorem ^ 
Theorem is proved in Section 5. Finally, we make some concluding remarks in 
Section 6. 



3. Upper Bounds on the Tail of the Wavepacket 

In this section we prove an upper bound on the tail of the (time-averaged) 
wavepacket in terms of transfer matrix norms for complex energies. It is shown 
how this leads to upper bounds on transport exponents in suitable cases. 

Observe that, by unitarity, J2 n a ( n ^) = 1 f° r am7 T. For an integer N > 1, 
define 

P r (N,T) = a(n,T) and P t (N,T) = ^ a{n,T) 

n>N n<-N 

so that P(N,T) = Pi(N,T) + P r (N,T). We call P r (N,T) (resp., P t (N,T)) the 
right (resp., left) outside probability. Our first goal in this section is to find upper 
bounds for P r (N, T) and P t (N, T). We will give the details for P r (N, T). Analogous 
bounds for Pi (N, T) can be obtained in the exact same way. 
To be precise, we want to find N(T) — > oo such that 

(29) lim P r (N(T),T) = 0. 

T — >oo 

The intuitive meaning of (|29|l is that for T large, the wavepacket ip(t) = e~ itH 5i, 
< t < T, has essentially no weight in the region (N(T), oo) for times up to T, 
at least in a time-averaged sense. A ballistic upper bound shows that (12911 always 
holds if we let N(T) = T 1+v for some v > 0. In the case of a dynamically localized 
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system (e.g., when the potential is given by independent identically distributed 
random variables), trivially holds for any N(T) such that N(T) — > oo. The 
most interesting case is the intermediate one, where (I29|) holds for N(T) — T a with 
some a £ (0, 1) (and does not hold for N(T) = T 7 with some < 7 < a). 

When proving (|29|l and the analogous bound for Pi(N(T),T), two kinds of results 
are of interest: 

(i) Find N(T) such that P r (N(T), T) and Pi(N(T),T) go to zero fast; for example, 
faster than any inverse power of T. In particular, if N(T) = T a with a G (0, 1), we 
get 

(30) fl+(p) <a+ <a, p>0. 

(ii) Find smaller values of N(T) such that P r (N(T),T) and P[(N(T),T) go to zero 
according to some inverse power of T. This allows one to prove upper bounds for 
a, and upper bounds for (3^ (p) that depend on p and are better than l|30|) . 

Here we are mainly interested in the first problem, which deals with the sharp 
front of the wavepacket. 

The following results gives upper bounds on outside probabilities in terms of 
quantities involving the norms of transfer matrices at complex energies. 

Theorem 7. Suppose H is given by Jjyl , where V is a bounded real-valued function, 
and K > 4 is such that o~(H) C [-K + 1,K— 1]. Then, the outside probabilities 
can be bounded from above in terms of transfer matrix norms as follows: 

P r (N,T)<exp(-cN)+T 3 J fmax^ ||$ (n, E + i) || 2 ^ dE, 

Pi(N,T) < exp(-cTV) + T 3 J ( ^max ^ ||$ (n, £ + ^) d£. 

T7ie implicit constants depend only on K . 

The starting point is the Parseval formula |24l 157) . 

1 /"°° 

(31) a(n,T) = — \{{H - E - ±,)-H Xl 8 n )\ 2 dE. 

Let us write e = 1/T and i?(z) = (H - zl)- 1 for zeC\l. We assume T > 1, so 
that < e < 1. Thus, we have 

e r 00 

(32) P r (7V, T) = — M r (N, E + ie) dE, 

* J -co 

where 

M r {N,z) = \{R{z)h,5n)? = IIxjv^^iII 2 , 

n>N 

and XAf( n ) =0: n < N, Xn{ti) = 1, n > JV + 1. The key problem is to bound 
M r (N,z) from above. If the energy E is outside the spectrum of the operator 
H, this is a trivial problem. Assume that E £ (—00,—K) U (K, +00), so that 
7] = dist(I2 + ie, cr(H)) > 1. In this case the well-known Combes- Thomas estimate 
yields 

(33) \(R(z)S 1 ,S n )\ < -exp(- C min{c77,l}|n- 1|) 

V 
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with some uniform positive c. Using (|33J) one can easily see that 

M r (N, E + ie) dE < C(K) exp(-dV), c> 0, 



' E:\E\>K 

and thus goes fast to for any N(T) — T a , a > 0. The main problem is to estimate 
L r (N,e)=[ M r (N,E + ie)dE, 

J-K 

where E + ie may be close to the spectrum of H. We will link M r (N,z) to the 
complex solutions to the stationary equation Hu = zu. 
Define two operators with truncated potential: 

H±i>{n) = ^(n - 1) + >4>{n + 1) + V* (n)V(n), 

where 

V^(n) = V(n), n < N, V^(n) = ±2K, n > N + 1. 

Such operators already appeared (in a different context) in [33]. We denote by 
R%(z) their resolvents and define 

7±fAT ~\ _ Ik. D± / „\ JT 1 1 2 



S ± (N,z) = \\ X nR%(z)6u 
Lemma 1. For any E S [— iT, K],0 < e < 1, we have 

e 2 S ± (N,E + ie) < M r (N,E + ie) < e^ 2 S ± (N, E + ie), 
where the implicit constants depend only on K . 
Proof. By the resolvent identity, it follows that 

R(z)S 1 - R%(z)S 1 = R(z)(V - V±)R±{z)8 x = R{z)( XN {V(n) T 2K))R±{z)8 x . 
Since ||i?(z)|| < e -1 , e < 1, and V is bounded, we get 

M(JV,z) < 2S ± (N,z) + 2\\R(z)(V-V±)R%(z)5 1 \\ 2 

< 2S ± {N, z) + 2e- 2 \\ XN (V T 2K)R%{z)5 1 \\ 2 

< c 2 (K)e- 2 S ± (N,z). 

The other bound can be proved in a similar way. □ 



The next step is to link the quantities 5 ( N, z) with the solutions to the sta- 
tionary equation. For any complex z, define uo(n, z) as a solution to Huq — 
zuq, lio(0, 2) = 0, uq(1, z) — 1 and u\(n, 2) as a solution to Hu\ = zu\, u\{Q,z) — 
f, u\{l,z) = 0. An important observation is the following: since V(n) — 
V^(n), n < N, the solutions uq,u\ for H and for coincide for all n < N + 1. 
We will consider uq, "i only for such n, thus we use the same notation for uq, u\ in 
the case of H and H^. 

For any E £ [-K, K], < e < I , z — E + ie, we define 

f 



= ^(zT2K+ ^(z T 2Ky-4). 

Observe that |z=F2i4T| > K > 4. Thus, one of the branches of the square root yields 
|Af(z)| < 1 and the other (z)\ > 1. For the values of z under consideration, we 
clearly have 

|A±(z)| < b(K) < 1, \X+(z) - \i(z)\ > c(K) > 
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with uniform constants b(K), c(K). The numbers X : 2 (^) in fact are the eigenvalues 
of the free equation 

(34) u(n-l)+u(n + l) = (zT2K)u(n), 
whose general solution is 

u{n) = CiAl 1 + C 2 A' 2 l . 

Lemma 2. Let z = E + is, where E e [-K, K], < e < 1. Then, for n > N > 3, 
we have 

(35) \(R%(z)8i,5 n )\ < 2e-i 

|Ai {z)u (N, z) - uq(N + 1, z)| 

and 

(36) | (fl ±(^ lA)| < £ -i ui 

|Aj (z)wi(iv, zj - ui(iV + 1, 2)| 

Proof. The following formula for the resolvent of any operator of the form holds 
(see, e.g., [ST]): 

(37) (R(z)5\, 5 n ) = d(z)uo(n, z) + b(z)ui(n, z), n > 1, 



(38) (R(z)Si, 8 n ) — d(z)uo(n, z) + c{z)u\(n, z), n < 1, 

where 

d(z) = 
b(z) = 
c(z) = 



— m + (z)m_ 


,{z) 


m + (z) + m. 


_{z) 


m_ (z) 




m + {z) + m. 


_(z) 


— m+(z) 





m+(z) + rti— [z) ' 

with some complex functions m+(z), m_(z) (which depend on the potential). Since 
11^(^)^111 < and 

d(z) = (R(z)8 1 ,8 1 ), c(z) = (R(z)S 1 ,8 ), 

we get \d(z)\ < e" 1 , \c(z)\ < e _1 . Since b(z) = 1 + c(z) and e < 1, we also get 
\b(z)\ <2e~ 1 . Summar lzmg. 

(39) \d(z)\<e-\\c(z)\<e-\\b(z)\<2e-\ 

One should stress that the bounds (|39|l hold for operators of the form (J3J) with any 
potential, and the constants are universal. 

Consider the operators and define (j) ± = Rjj(z)8i (the vectors <jr^ depend 
on N, of course). Since 

the function (^(n) = {(^,8^) obeys the equation H^4>(n) = z4>(n),n > 2. Since 
V]y(n) = ±2K, n > N + l, cj>(n) obeys the free equation for n > N+ 1. Hence, 

(40) (cb ± (N + k + l),<b ± (N + k)f = C 1 ± A±(z) fe ei + C 2 ± A 2 t (z) fe e 2 , fc > 0. 
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Here ef 2 — {-^i2( z )t^) T are two eigenvectors of the matrix corresponding to the 
equation (|34l) , and the constants , are defined by 

(^(N+l), ^(AOf = Cfex + C±e 2 . 

Since \Xf (z)\ < 1, |A^(z)| > 1 and ± e ^ 2 (Z), the identity igOl implies that 
= 0, and thus 

(41) (^(N + l^^WF = D%(z)(Xf(z),lf, D±(z)^0. 
On the other hand, 13 7() implies 

(42) ^(N + 1) = d%(z)u (N + 1, z) + &£(«)ui(JV + 1, z), 



(43) ± (iV)=4(2)« o (iV,z) + 6^(^i(iV^) ) 

Remark that d^(z) ^ 0, since it is a Borel transform of the spectral measure corre- 
sponding to the vector <5i [5 7) . The same is true for b N (z), c N (z) since m + (z), m_ (z) 
are functions with positive imaginary part (it also follows from d^(z) ^ and the 
expressions of b, c) . Using the fact that uo(n + l, z)ui(n, z) — uq{ti, z)u\(n + l, z) = 1 
for any n (in particular, for n — N), and (|41() - (|43|l . it is easy to calculate Dj^(z): 

(44) Dt(z) = ^ = ^ 

V ; NK ' Xf{z)ui(N,z)-Ui(N + l,z) Xf(z)u (N,z) - u {N + l,z)' 

As observed above, the solutions uo, u i are the same for H,H^ if n < N + 1. It 
follows from gUJ), where Ci = D ± and C 2 = 0, that 

(R±(z)6 1 ,6 n ) = D±(z)Xf(z) n - N , n>N. 

The result of the lemma follows now directly from 1)39(1 and (|44|) . □ 

Lemma 3. for z = E + ie with E 6 [—A", if] and < e < 1, and N > 3, we /iai/e 

-l 

2 ' 



M r (N, z) < C(K )£" 4 max ||$(n,z) 

\3<n<iV 

Proof. The second bound of Lemma H an d the bound Ij35|) of Lemma [5] yield 

oo 

(45) M r (N, z) < A(K)e- 4 |A±(z) Wo (^ z) - u (N + 1, z)\~ 2 £ |A±(z)| 2fe 

(46) < B(K)e^\Xf(z)u (N,z)-u (N+l,z)\- 2 

with uniform B(K), since |Af (z)| < b{K) < 1. 

Denote E/* = Xf(z)uo(N, z) — uq(N + 1, z). It is not hard to see that 

(47) \U-\ + \U + \>\\- l {z)-\i{z)\\u {N,z)\ 
and 

(48) \U-\ + \U+\ > \\T(z)\\U-\ + \Xt(z)\\U+\ > \Xi(z) X+(z)\\u (N + l,z)\. 
Since \X±(z) - X+(z)\ > c(K) > 0, (g^-lgHJ imply 

M r (iV, z) < C(if) £ - 4 (|w (iV, z)\ 2 + \u (N + 1, z)^)' 1 . 
Using l(36|l . we can prove a similar bound with uq replaced by u%. Thus, 

M r (N,z) <C(K)e- 4 \mN,z)\\- 2 . 
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Since M r (n, z) is decreasing in n, the asserted bound follows. □ 

Proof of Theorem^ The assertion is an immediate consequence of (|32|) and 
Lemma|31(and the analogous results on the left half- line). Note that we can replace 
[3, N] by [1, N] since this modification only changes the ^-dependent constant. □ 

Proof of Theorem^ Let us choose N(T) = [CT a \, where C € (0,oo) and 
a e (0,1) are chosen such that JT3J and O hold. Observe that P(N(T),T) = 
P([N(T)\,T). Then Theorem □ shows that P r (N(T),T) and Pi(N(T),T) go to 
faster than any inverse power of T. By definition of S + (a) and cn+ (cf. Q and 
©), it follows that a+ < a, which is JTSJ. Finally, JTJjJl follows from ptljl. □ 

Let us now show that the transport exponents a„ can be expressed in terms of 
the integrals 

I(N,T) = J (j\<f>(N,E + ±)\\~ 2 + ||$(-iV,£ + i)||" 2 ) dE. 
It follows from Theorem [3 that 

(49) P(N, T) < exp(-ciV) + T 3 I{N, T). 

Define 

log/(LT Q -2j,T) 



W (a) = — lim inf 
W + {a) = — lim sup 



T^oo logT 

log/(LT°-2j,T) 
/ v logT 



7 ± = sup{a > : W ± (a) < oo}. 
Theorem 8. We have that 

a± = 7 ± 

Proof. We first prove a lower bound for a(n + 1,T) + a{n,T). Let = 
{R(z)5i,8 n ), where z = £ + i/T. The identities j57^. imply 

(0(n + 1), 0(n)) T = $(n, z)(d(z), 6(^)) T , n > 1, 
(0(n + 1), 0H) T = *(n, z)(d(z), c(z)) T , n < 1. 

Since det$(n, = 1, we obtain (see, e.g., [23] for details) 

(50) \<j>{n + 1)| 2 + \4>(n)\ 2 > \mn,z)\\- 2 \d(z)\ 2 

for any n. Since d(z) — (R(z)Si, Si) is the Borcl transform of a compactly supported 
measure, it is easy to see that \d(E + y)| > lmd(E + y) > C/T with a constant 
that is uniform in E e [-K,K],T> 1. Therefore, J3TJl and (JSDJl yield 

C f K 

a(n+l,T)+a(n,T)> —g / ||$ (n, £ + 4,) |f 2 dE, 



and 
(51) 



P{N,T)>-^J (||$ (TV + !,£+!) || 2 + ||$(-iY-2,£+i)|| </£. 
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Since the potential V is bounded, and E £ [—K,K], T > 1, we have that 
\\$(N + l,E+i)\\ x \mN,E+±)\\, 
\\*(-N-2,E+i)\\ x \m-N,E+±)\\, 
with constants depending only on K. Therefore, l|51|) implies 

(52) P{N,T)>^I{N,T). 

Let us show that a+ = 7 + . Assume first that 7 + < a+. Then there exists a such 
that 7 + < a < a+. The definition of 7 + implies PF + (a) = oo and thus 

(53) I([T a ~2\,T)<^. 

for any M > uniformly in T. It follows from igHJ that P(T Q - 2,T) = 
P([T a - 2j,T) also obeys a bound like and thus 5+ (a) = oo. However, 

this is impossible since a < a+. Therefore, 7 + > a+. Similarly, 7 + < a J follows 
from the bound 152|) . The identity a~ = 7~ can be proven in the same way (with 
inequalities like (|53H valid for some sequence of times Tk — > 00). □ 

Proof of Theorem^ It follows directly from (|17fl and the definition of W ± (a) that 
IWfta) - W^WI < 2 A for any a, and hence 7^ =7^. The result now follows 
from Theorem^ (Since < 1, it is sufficient that i|17|) holds for \n\ < 1/e.) □ 

4. The Fibonacci Operator 

In this section we prove upper bounds for the transport exponents f3 + (p) in the 
case of the Fibonacci potential given by (|18fl for A sufficiently large. 

In order to apply Theorem \7\ we need to find lower bounds for the norms of 
transfer matrices at complex energies. Rather than norms, we will study traces of 
transfer matrices. This will be sufficient because lower bounds for traces give rise 
to lower bounds for norms in a trivial way. Moreover, the study of traces is natural 
in the case of the Fibonacci model since it displays a well-known hierarchical struc- 
ture on the level of transfer matrices that induces recursion relations for the traces 
of transfer matrices over intervals whose length is given by a Fibonacci number. 
This in turn is most conveniently described by a dynamical system, the so-called 
trace map. We remark that in the existing proofs of lower bounds for transport 
exponents, it was required to find upper bounds for transfer matrix norms. Inter- 
estingly enough, it was sufficient also in that case to study traces. This less trivial 
transition from traces to norms can be achieved in the situation at hand using an 
observation of Iochum and Testard 07]; see also |46|. 

While the trace map has been studied by many authors (e.g., [TT1 IT71 1571 1551 
15^1 IH51 EH), it has been analyzed so far only for real energies. Thus, the trace 
map dynamical system is usually considered as a map from R 3 to M 3 . Theorem UJ 
suggests studying complex energies. Thus, we are led to study the trace map as a 
complex dynamical system, that is, as a map from C 3 to C 3 . The initial element of 
C 3 will be an energy-dependent vector. We investigate the stable set, that is, the 
set of complex energies such that the associated vector in C 3 has a bounded trace 
map orbit. This set of energies is, in fact, a subset of K, but it is important to 
study the canonical approximants to this set as subsets of C. Their complements 
are escaping regions in the sense that once an orbit enters such a set, it escapes to 
infinity at a super-exponential rate. Since all non-real energies give rise to escaping 
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trace map orbits, we can study in this way the maximum number of iterates it 
takes, for a given imaginary part e of the energy z, to enter an escaping region. 
From this point on, we can control the increase of the trace and this will prove 
sufficient for us to obtain the lower bounds we need. 

Our goal is to derive Theorem [3] from Theorem ^ Since V{— n) = V(n — 1) for 
n > 2 (see |fi9|). we can restrict our attention to one half line. Thus, we will give 
details only for with suitable a. The proof of (|14|l is completely analogous. 

For z £ C, define the matrices M k (z) by 

(54) *{F k ,z)=M k (z), k>l, 

where F k is the fc-th Fibonacci number, that is, Fo = F% = 1 and F k+ i = F k +F k _i 
for k > 1. It is well-known that 

Mfc+ifa) = M k -x{z)M k {z), k>2. 

For the variables x k (z) — (Tr M k (z))/2, k > 1, we have the recursion 

(55) x k+ i(z) = 2x k (z)x k -i(z) - x k -2{z) 
and the invariant 

A 2 

(56) x k+ i(z) 2 + x k (z) 2 + x k -i(z) 2 - 2x k+ \(z)x k (z)x k -i(z) - 1 = — . 

Letting x-i(z) = 1 and xq(z) — z/2, the recursion 1)55(1 holds for all k > 0. See, for 
example, [^1^12^1) f° r these results. 
For 6 > 0, consider the sets 

o{ = {z£ C : \x k {z)\ < 1 + 5}. 

For S C C, we let S R = S n E. The set erg is the spectrum of the fc-th periodic 
approximant to the Fibonacci operator. In other words, x k is the discriminant of 
this periodic operator and hence the set (<J k ) R consist of F k closed bands that, 
in general, have disjoint interiors. However, Raymond has shown that in these 
particular cases, even the closed bands are mutually disjoint |65) . Thus, there are 
exactly F k disjoint bands making up the set (<7 k ) R , and each one of them contains 
exactly one zero of x k and each one of the F k — 1 bounded gaps contains exactly 
one critical point of x k . 
We have 

4 u 4-i ^ 4+i u 4 -> 

the latter set being the spectrum of the Fibonacci Hamiltonian. Outside of these 
sets, the traces grow super-exponentially. More precisely, a modification of Siito's 
proof shows the following: 

Lemma 4. A necessary and sufficient condition that x k (z) be unbounded is that 

(57) \x N -i(z)\ < 1 + 6, \x N (z)\ > 1 + 6, \x N+1 (z)\ > 1 + 6 
for some N > 0. This N is unique. Moreover, in this case we have 

(58) |ar n+2 (z)| > \x n+ i(z)x n (z)\ for n> N 
and 



(59) 



\x N +n(z)\ > {1 + S) F - forn>0. 
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FIGURE 1. The structure of the sets (cr| ) R when (J^TJ holds. 

Proof. Suppose that l|57[) holds true for some N > 0. Then 

\xn+2(z)\ > \x N +i(z)x N {z)\ + (\x N+ i(z)x N (z)\ - \x N -i(z)\) 

> \x N+1 (z)x N (z)\ + 8(1 + 8) 

This also implies \xn+2(z)xn+i(z)\ > \xn(z)\. Thus, 1(58(1 follows by induction. 
This in turn shows both 1(59(1 and that there is at most one N > with (|57|l . 

Now suppose that (|57|l holds for no value of N > 0. Since x_i(z) = 1, this 
implies that for every n with |a;„(z)| > 1 + 8, we must have both |x„_i(z)| < 1 + 8 
and |x„+i(z)| < 1 + 8. The invariant i(56(l therefore shows that the sequence x n (z) 
must be bounded. □ 

Moreover, assuming A > Ao((5), where 

(60) X (S) = [12(1 + <5) 2 + 8(1 + 8) 3 + 4] 1 / 2 , 
the invariant 1)56(1 implies 

(61) 4 n 4+i n 4+2 = 0- 

Lemma 5. Fix 8 > and A > A (<5). Then the set (t^) r consists of closed 
bands that are mutually disjoint. Each one of them contains exactly one zero of X}. 
and each one of the — 1 bounded gaps contains exactly one critical point of Xk ■ 

Proof. For <5 = 0, these statements are known and were recalled above. If <5 > 0, let 
t run from to 8. The bands of (cr£) M fatten up and what needs to be avoided is 
that two consecutive bands touch. By the assumption A > Ao(<5) this is impossible, 
however, as can be seen from (|61|l and the structure of the sets (4 ) H induced by this 
property (compare Figure 1 and the discussion of the sets (o^) R in 57. Section 5] 
and H21 Section 6]). □ 

Lemma 6. Fix 8 > and A > Ao(<5). Then the set 4 nas exactly Fk connected 
components. Each of them is a topological disk that is symmetric about the real 
axis. 

Proof. Since Xk has degree Fk, the set 4 ^ as at most Fk connected components. 
Since the coefficients of Xk are real, it follows also that each such component must 
be symmetric about the real axis. On the other hand, Lemma El shows that the set 
{4 ) K h as exactly Fk connected components. Thus, the maximum modulus principle 
show that each component of 4 lii a topological disk and there are exactly Fk of 
them. □ 
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The final ingredients we need are the following bounds on |a/ fc (z)| for z £ (c^) R 
(cf. (S3 Proposition 5.2] and Equation (57)]): 

Proposition 1. If A > 8, then 

(62) w k (z)\>m k/2 

for all fc > 3 and z £ (cr°) R . Here, 



^ (A) = A-4+ V(A-4)2-12 

and hence £(A) = A + 0(1) as A — * oo. 
Proposition 2. If X > 4, i/ien 

(63) < C(2A + 22) fc 
/or aH fc > 1 and z G (cr°) R . 

Now we can prove our main result on the structure of the set af.. Denote 
B(z, r) — {w £ C : \w — z\ < r}. 

Proposition 3. Fix k > 3, <5 > 0, and A > Xq(2S). Then, there are constants 
c$,d$ > such that 

(64) [jB(zi j \r k ) £a s k C Qs(z«,i? fc ), 
i=i i=i 

w/iere {z£ }i<j<F k are the zeros of xu, r k — c$(2X + 22)~ fe j and Rk = dgti(X)~ k t 2 . 

Proof. Consider a connected component Cj of a\ s . By the assumption A > Aq(25), 
Cj contains exactly one of the Ffc real zeros, z^\ Moreover, it contains exactly one 
connected component of af. , which we denote by Cj . Our goal is to show that with 
the radii rfc , Rk described in the proposition, 

(65) B (z^ , r fc ) C Cj C S (z<f } , ife) . 

Clearly, (|51J follows from 

Since Cj contains exactly one zero of Xk, it follows from the Maximum Modulus 
Theorem and Rouche's Theorem (e.g., Theorems 6.13 and 10.10]) that 

x k :int(Cj)-> B(0, 1 + 2<5) 

is univalent, and hence 

x^ 1 : B(0, 1 + 25) mt(Cj) 

is well-defined and univalent as well. Consequently, the following mapping is a 
Schlicht function: 

, N , , x7 1 ((l + 25)z)- z{ j) 

F:B(0,1)^C, F(z) = k vv — '~4 — • 

{ 1 {) (l + 2o)[(x,T 1 )'(0)] 

That is, F is a univalent function on B(0, 1) with F(0) = and F'(0) = 1. 
The Koebe Distortion Theorem (see ^3 Theorem 7.9]) implies that 

^ JTTW ~ lF{z)l ~ (T^W f ° r N " L 
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Evaluate the bound l|66|) on the circle \z\ = jr^g- For such z, we obtain 

jl + 5)(l + 25) (l + 5){l + 25) 

(2 + 35)2 S|* WIS p 

By definition of F this means that 

\ X - k \{l + 25)z) - *«| < {1 + S) £ + 2S \ l + 25)|K 1 )'(0)| 

and 

I^Hd + 26)z) - ,«| > + 2<5)|(x- 1 )'(0)| 

for all z with |z| = ^r^. In other words, if |z| = 1 + 6, then 

(67) I^W-^l< (1+ ^ + 2i)2 |(O'(0)l 
and 

(68) K i W -,W|>(^^, (jB -i m | 

Since |(a;^ 1 )'(0)| = |a;' fc (z^)| _1 , we obtain from Proposition ^ and (joTjl 

(69) K 1 W-4°l<( (1 + ^ 1 + M) )%(A)- fc/a 
for all z of magnitude 1 + 5. Similarly, Proposition and Ij68(l give 

for these values of z. Note that as z runs through the circle of radius 1 + 5 around 
zero, the point (z) runs through the entire boundary of Cj. Thus, (|65|) follows 
from and fijfo. □ 

Proof of Theorem^ By the assumption A > 8, it is possible to choose 5 > such 
that A > Xq(25); compare ||dT7|) . Since the z^ are real, the last proposition gives in 
particular that 

<? S k C {z e C : |Im z| < rf<5€(A)- fc/2 } C {z e C : |Im z| < d a F fc ~ 7(A) }. 
for a suitable 7(A). This implies that 

(71) 4 U 4 +1 C {z e C : |Im z\ < d s F^ W }. 

If r] = + l)/2 (so that Fk behaves like r/ k for k large enough), then the 
constant 7(A) may be chosen according to 

= logg(A) 
71 J 2(1 + 1/) log 7? 

for some ^ > 0. Since we are interested in large T behavior, we can work with any 
positive v and then have the inclusion above for k > ko(v). 

In other words, for each e = Im z > 0, one obtains lower bounds on \x n (E + ie)\ 
which are uniform for E G [— JT, K] C R. Namely, given e > 0, choose fc minimal 
with the property dsF k ~ J ^ < e. By ifTTJ) . we infer that |a:jfe(JS + ie)| > 1 + 6 and 
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\xk+i{E + ie)\ > 1 + 6. Since \x~i(E + ie) \ = 1 < 1 + 5, we must have the situation 
of Lemma 01 for some N < k. In particular, for n > k, H59JI shows that 

\x n {E + ie)\ > {l + 6) F - k . 

This motivates the following definitions. Fix some small 6 > 0. For T > 1, 
denote by k(T) the unique integer with 

P 7(A) jpfW 
r k(T)-l < rp r k(T) 



and let 

N(T) = ^fc(T)+[ % /fc(T)J • 

Thus, for every v > 0, there is a constant C„ > such that 

(72) iV(T) < CoT^T . 

It follows from Theorem and the argument above that 

PJN(T),T) <exp(-cN(T)) +T 3 [ ( max 11$ (n,0;E + i) ll 2 ] dE 

J-K \3<n<N(T) 1/11 J 

<exp(-c7V(T))+T 3 (l + (5) lVm^j. 

From this bound, we see that P r (N(T),T) goes to zero faster than any inverse 
power. We note again that we can obtain the same result on the left half- line due 
to the symmetry of the potential. Therefore, by (|72[) . we can apply Theoremnjwith 

1 2(l + ^)log» 

7(A) log£(A) 

and hence (|23[1 and (|24|) follow from l|15|l and (|16f) . respectively, since we can take 
^ and £ arbitrarily small. □ 

5. The Almost Mathieu Operator 

Suppose we are given an operator family {Hg^} satisfying the assumptions of 
Theorem |SJ 

Let K = H/Hoo + 3. Then, 

(73) (r(H e<u ) C [-K+1,K- 1] for all u € [0, 1). 
We will show the following 

Proposition 4. Assume that Hg u obeys the assumptions of Theorem\^ Denote 

I r (N,T)= ( C max ||$ (n, E + i, 0, w) || 2 ) dE, 

Ii(N,T) = J ^max^ ||$ (rc, E + j,, 0, uj)\\ 2 ^j dE. 

(a) For every a > 0, i/iere exists a sequence Tk — * oo suc/i i/iai 

(74) 7 r (T^T fe ) = 0(T fc —) 
and 

(75) I l (T£,T k ) = 0(T^) 
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for every m > 1. 

(b) Assume moreover that 8 obeys Ij27(l . Then for any a > 0, 

(76) I r (T a ,T) = 0(T- m ), Ii(T a ,T) = 0(T~ m ) 
for every m > 1. 

Remarks, (a) The a-dependent sequence {Tk} will be related to the denominators 
qk of the continued fraction approximants Pk/qk of 9 (cf. |56|L 
(b) Our proof of Proposition 0] will have some elements in common with the argu- 
ments given by Jitomirskaya and Last in |52| in their proof of zero-dimensionality 
of spectral measures. 

Proof of Theorem^ The assertion is an immediate consequence of Proposition 0] 
together with Theorem ^ and the first remark after that theorem. □ 

Proof of Proposition^ Denote 

K = {z £ C : -K < Rcz < K, < Imz < 1}. 
It is easy to see that there exists a constant T' such that 

(77) \\®(n,z,6,u))\\ < e |n|r ' for all u £ [0,1), n£Z, z£K. 
Combining this with the assumption of Theorem we obtain 

(78) < T < 7e (z) < T' < oo for all z £ K. 
For n > 1, let 

'n-fe{zY 



A n = Uo£ [0,1) : \\§{n,z,6,u)\\ > exp 
It follows from and J77|) that 

nie(z) < / log||$(n, z, 9,u)\\ du 



]og\\$(n,z,9,u)\\dw+ / log||$(n,*,0,w)||dw 

A n J{0,l)\A n 

< \A n \nT' + (l-\A n \)n^-, 

where | • | denotes Lebesgue measure. 
Therefore, 

1 ~ 2T'- 79 (z) ~ 2r-T~ 
Let d be the degree of the trigonometric polynomial f. 8 Then each entry of 
<I>(n, z, 6, w)*$(n, z, 0, w) is a trigonometric polynomial of degree at most 2nd, and 
the set A n consists of no more than And intervals. Therefore, there exists an interval 
/„ C A n with 9 

w >- 53- 

Let 

n * = Lid J + 



^Note that we have d > 1 since the assumptions of the theorem preclude the case d = 0. 
9 If one considers the Hilbert-Schmidt norm 1 1 <E» 1 1 ^ = Tr(<3?*<3?), which is of course equivalent to 
the operator norm, this claim is obvious. 
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By a standard lemma in continued fraction approximation (see, e.g., |52l Lemma 9] 
for a formulation suitable to our purpose), for every ui, there exists a j £ 
{0, 1, . . . , qk + qk-i — 1} such that j6 + oj mod 1 belongs to I Uk , and hence to 



A nk . That is, ||*(n fc ,M,.70 + w)|| >exp(^f^). 



Since 



$(j + n k ,z,6,u;) = &(nk,z,6,j0 + w)$(j,z,0,w) 
and each $ is unimodular, we see that either + n h , z, 0, or z,6,uj)\\ 



is greater than exp 
Let 



n k ~fe(z) 
4 



jk = min { j G {0, . . . , qk + Qk-i - 1 + n k } : ||$(j, z, 0, w)\\ > exp 



V 4 

By definition of jk, (|77|l . and l|78|) . we have < jk and hence there are z- 
independent constants C\, C 2 such that 

(79) C iqk < j k < C 2 q k . 
Moreover, there is a z-independent constant C3 such that 

(80) \\$(j k ,z,6,uj)\\>exp(C 3 qk) 
Thus, the definition of I r {N,T) and imply 

(81) I r (C 2 qk,T) < 2KeM-2C 3 q k ) 

for any T > lieN. To prove part (a) of the proposition, consider, for a given 
a > 0, the sequence T fc = [C 2 q k ] 1/a . The bound {ZU follows from JSTJ. Since 

(82) $(n,z,(9,w) = $(-n,z,-0,0 + w) 

and the q k s are the same for 9 and — 0, we see that we can prove (|75|l with the 
same sequence {T k }. 

Let us now prove part (b) of the proposition. Choose some a > 0. For any given 
T > 1, one can find some k = k(T) such that 

(83) (C 2% ) 1/Q < T < (C 2 qk+i) 1/a - 
Let 

(84) 7 fe = . 

qk 

It follows from the definition of 7 r , (jSTJ) and ||S3}-||HI|! that 

I r (T a ,T) < I r (C 2 q k ,T) < 2Kcxp(-2C 3 q k ) = 2Kq k l^ hk < C T- 2aC3 ^. 

By assumption we have lim^oo 7^ = 0. Since fc(T) — > 00 as T — > 00, we get 
I r (T a ,T) = 0(T- m ) for every m > 1. Due to the observation (JH2J, the same is 
true for Ii(T a ,T). This gives (|76|l and concludes the proof. □ 
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6. Concluding Remarks 

We introduced a general upper bound on outside probabilities in terms of trans- 
fer matrix norms at complex energies. We remark that a completely analogous 
proof shows the same result for operators of the form ((3J) in £ 2 (Z + ) with any self- 
adjoint boundary condition at zero (e.g., Dirichlct). This upper bound for outside 
probabilities can be used to bound the transport exponents from above in certain 
situations. 

While we carried out the analysis for the initial state 8\ , one can prove dynamical 
upper bounds for more general initial states. Namely, one can take ip = f(H)8i, 
where / is some smooth compactly supported function, and apply ideas from |72j . 
In this way one can study dynamics locally in the energy. For example, if H has 
some absolutely continuous spectrum, but the spectrum is purely singular on [a, b] , 
one can consider functions / that are supported on [a, b], so that ip belongs to the 
singular subspace of H . 

The general result was used to prove quantum dynamical upper bounds for the 
Fibonacci operator in the large-coupling regime. This constitutes the first example 
for which a non-trivial upper bound on flf (p), < p < oo, can be established 
outside the realm of dynamically localized systems. In particular, we obtain the first 
rigorous anomalous transport result. The asymptotic dependence on the coupling 
parameter follows a law that was numerically predicted by Abe and Hiramoto in 
1987. 

Combining the ideas of the present paper with those from it is possible 
to perform a similar quantum dynamical analysis for general Sturmian potentials 
given by 

Moreover, it is possible to improve the known dynamical lower bounds for the 
Fibonacci Hamiltonian by employing the first inclusion in (|64[l together with the 
main result of |24| . 

We also showed that all lower transport exponents and, under a weak assump- 
tion on the frequency 9, all upper transport exponents vanish for quasi-periodic 
potentials in the regime of a uniformly positive Lyapunov exponent, that is, for 
sufficiently large coupling by a result of Herman. For the particular case of the 
almost Mathieu operator, we obtain this result in the regime A > 2. 

We expect that our method will be useful also in the context of the almost 
Mathieu operator at critical coupling, A = 2, where almost-diffusive transport is 
expected. It is an interesting open problem to bound the transport exponents in 
this case. For Diophantine frequencies, Bellissard, Guarneri, and Schulz-Baldes 
have established a lower bound on phase-averaged tranport exponents in terms of 
the multifractal dimensions of the density of states ^Oj- Since there are no known 
lower bounds for these dimensions, however, this result does not yield any explicit 
estimates for the dynamical quantities at this point in time. 
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